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Abstract 

We construct a non Abelian model for SU(2) QCD in Euclidean 
three-dimensional spacetime and study its different phases. The model 
contains a center vortex sector coupled to a dual effective field encoding 
information about how the vortices are paired in the ensemble. The 
possible phases in parameter space are interpreted in terms of the 
proliferation of either closed center vortices or closed chains, where 
the endpoints of open vortices are attached in pairs to monopole-like 
defects. 

1 Introduction 

The vortex model introduced by t' Hooft |T] is a low energy effective theory 
that successfully describes some aspects of the confinement mechanism in 

2 + 1 dimensional SU(N) Yang-Mills theories. It is defined in terms of a 
dynamical variable which is a complex scalar field V equipped with a discrete 
Z(N) symmetry, realized with a Lagrangian, 

C = d^Vd/jV + [i 2 VV + aiyVf + /3(V N + V N ) . (1) 

Its form is based on a study of the possible nontrivial vortex correlation 
functions in the original theory. In particular, the confining phase is described 
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as one where the discrete Z(N) symmetry is spontaneously broken, due to 
the presence of a vortex condensate. The possibility of an effective field 
representation for 3D center vortices relies on the fact that an ensemble 
of stringlike objects can be thought of as a sum over different numbers of 
particle worldlines, which corresponds to a second quantized field theory [2]- 
[JJ. Based on this observation, one of us proposed, in refs. [HIE], a generalized 
vortex model where <9 M in eq. (jl]) is substituted by the covariant derivative 
Dp, that depends on a dual vector field A M describing the off-diagonal sector. 
The dynamics is completed with a Proca action term for A M . In ref. [10J, we 
derived this model by considering an ensemble of chains, where the vortex 
endpoints are attached in pairs to monopole-like defects, and following recent 
polymer techniques to compute the vortex end-to-end probability. 

Effective field models can also be obtained in scenarios based just on the 
monopole (instanton) component. In this case, the assumption of Abelian 
dominance and the associated monopole ensemble is encoded in a sine-Gordon 
type model for a scalar dual field (see fl2] and references therein), as oc- 
curs in the case of compact QED(3), discussed by Polyakov in ref. [13J. 

In spite of the fact that the initial theory is a non Abelian one, these 
effective models are Abelian, that is, additional information regarding this 
transition is already incorporated, while it would be desirable to see it ap- 
pearing as a phase transition in a previous non Abelian model. 

In addition, the different ideas regarding the magnetic sector have been 
explored in the lattice, relying only on monopoles |14|-|17|. only on center 
vortices |18|-|22j. or on chains [23j-[25j. Therefore, it would be interesting to 
construct a model where the possible phases in parameter space correspond 
to the different ensembles. 

In this article, we construct a non Abelian effective model which en- 
compasses a description of interacting effective gluons and center vortices. 
Depending on the choice of parameters, the vortices can be found in differ- 
ent states, including a phase where they are closed, and a phase where their 
endpoints become paired to form closed chains. 

To that aim, we use a parametrization that treats the different color 
components in a symmetric way [26J, and describes correlated monopoles 
and center vortices as defects of a local color frame h a , a = 1,2, 3. This 
parametrization is based on the usual manner to introduce thin center vor- 
tices in Yang-Mills theories [271 [28] . and corresponds to a symmetric form 
of the Cho-Faddeev-Niemi (CFN) decomposition |29|-|31|. used to represent 
monopoles as defects of the third component n = 77,3. For a description 
of center vortices in the CFN framework, and related consequences in the 
continuum, see refs. 0,12] • 

Since center vortices can be joined in pairs to pointlike monopoles, the 
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natural non Abelian field content of the model is given by a scalar field 
with one (magnetic) color index, generalizing the vortex field V in eq. ([!]), 
and a scalar field with two color indices, generalizing the scalar dual field in 
scenarios only involving the monopole component Q The order parameters 
present in the effective model bear a relation to the nature of the phase 
transition one may describe. In this respect, the interesting point has been 
raised |12| about whether the confining/deconfining phase transition is of 
the KT or Ising model type. The former involves the monopole sector: at 
high temperatures, the instanton magnetic flux is distributed along the two 
spatial directions, thus leading to effective logarithmic interactions. Then, 
because of dimensional reduction, instantons and anti-instantons tend to be 
suppressed by forming pairs. On the other hand, the latter naturally involves 
the vortex degrees of freedom, as they are the objects where the discrete 
symmetry transformations act. 

In the model we construct and study below, since it does contain order 
parameters for both the center vortices and the distribution of monopole- 
like defects they can concatenate to form chains, an interesting framework 
to discuss the competition between different phases shall emerge, originating 
a phase diagram with a rich structure. This paper is organized as follows: 
in section |2j we deal with the topological defects included in the model, in 
particular, their parametrization, and the functional and ensemble integra- 
tion over them. In section |3j based on the previous section results, and after 
discussing the possible symmetries, we construct an action for the effective 
model in terms of the fields introduced therein. Finally, in section |4j we 
present a study of the phase structure of the model, based on some assump- 
tions about the relative strength of its different terms. 



2 Non Abelian defects in YM theories 

We shall start from the SU(2) Yang-Mills action, Sym, which may be written 
as follows: 

Sym = \Jd 3 x F„ v ■ , (2) 

where F^ v is the non Abelian field-strength tensor. We use an arrow on 
top of any object to denote the 3-component vector formed by its compo- 
nents on the su(2) Lie algebra basis, whose elements are the (Hermitian) 
generators (T a )^ =1 . In the concrete case we are considering, they can be 
conveniently realized as T a = r a /2, where r a denotes a Pauli matrix; they 
satisfy [T a } T b ] = ie abc T c , and tr(T a T b ) = \5 ab . 

1 Interestingly, isospin two order parameters appear in models for liquid crystals |32| . 
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Thus, with this notation, we may write down the defining equation for 
Fpv, as follows: 

F^-f = - [D^D V ) , Dp — dfj, — igA^ ■ f , (3) 

where has been used to denote the covariant derivative operator, when 
acting on fields in the fundamental representation. 

It goes without saying that a 'canonical color basis' [e a ) 1 (with color 

components = 5 a &) can be introduced, so that = A ^ e a . This seemingly 
trivial remark is made in order to highlight the next step; namely, that one 
could have used a different basis. Indeed, in order to describe configurations 
with defects, in a symmetric way that admits its extension to finite size 
objects, we introduce a space-dependent color basis (n a )^ =1 , related to the 
original one by: S'T a S' _1 = n a ■ T (with S G SU(2)). Thus the new basis 
is connected to the canonical one by an orthogonal space-dependent matrix 
R(S): h a = R(S)e a , which belongs to the adjoint representation. In this 
representation, the corresponding infinitesimal generators shall be denoted by 
M a , with (M a ) bc = -ie abc . They satisfy [M a , M b ] = ie abc M c , tr (M a M b ) = 
25 ab . At this point, and equipped with the local basis, we consider the 
parametrization of the gauge field [26J: 

K = (a; - c;) n a , (4) 

where the frame dependent fields, 

C; = -j- g e abc h b .d,n c , (5) 

satisfy the properties: 

n b ■ 8^n c = -ge abc C; , C«M a = l -R- l d^R. (6) 

This corresponds to a symmetric form of the Cho-Faddeev-Niemi (CFN) 
decomposition [29j-[31j. In terms of the parametrization (|4]) of the gauge 
field, we note that the field-strength tensor becomes: 

Fft V = G* u h a , G° v = F* V {A) - J^„{C), (7) 

with F^A) = d^Al — d v A a ^ + ge. abc A b ^A c u (and an analogous expression for 
J~® U (C)), while the Yang-Mills action is given by: 

Sym = J d 3 x \g; u G; u . (8) 
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Regarding the color components of the frame-dependent tensor J~u V (C), they 
can also be obtained by commuting covariant derivatives in the adjoint rep- 
resentation: 

J*,(C) M a = - [V„ V v \ , V^d,- igC;M a , (9) 
so that the second equality in ^ implies the alternative expression for 

•?>(C) = y 9 tr i^R' 1 ^ d M ■ ( 10 ) 

This equation highlights the meaning of J-"^,(C), by showing that it can only 
be different from zero where R has defects; these, are characterized here 
by the noncommutativity of the mixed partial derivatives. These defects 
are zero measure objects; in other words, the partial derivatives will fail 
to commute on zero measure regions. Being this an effective theory, this 
should be interpreted as the assumption that the model describes physics at 
distances much larger than the size of the defects. 

Of course, there are infinitely many different local frames, and corre- 
sponding fields A a ^ that can be used to describe one and the same gauge 
field configuration, A a . One can use that large amount of freedom in order 
to split it into its 'regular' and 'singular' parts. Indeed, the A't measure 
will represent topologically trivial fluctuations. The singular configurations, 
described by the frames, will have a measure representing an ensemble inte- 
gration over defects. 

In ref. [2E], one of us has shown that the configuration in Q is tantamount 
to the usual way fZT\ |2"5] to introduce thin center vortices on top of a trivial 
field configuration «4^e a , namely, 

Ay, ■ T — SAfj, ■ TS- 1 + -Sd^S' 1 - I^S) ■ f . (11) 

Because of the presence of the last term, this is not just a gauge trans- 
formation of the topologically trivial gauge field. Indeed, the I^(S) field 
corresponds to the so called ideal center vortex, and is localized on a hyper- 
surface S. This is the region which, when traversed, makes S change by a 
center element. It is designed to cancel the contribution in the second term 
originated from the discontinuity of if? -1 , only leaving the effect of the border 
of S where the thin center vortices are located. That is, we can write, 

-sd ll s- 1 y = i ll (s)-f, (12) 



5 



where the subscript in the left-hand side amounts to just keeping in the 
calculation the term originated from the derivative of the discontinuity in 
S^ 1 . Considering two regular mappings U, U, the ideal vortex satisfies, 

I^USU- 1 ) ■ f = UI^S) ■ TU~\ (13) 

obtained from IJi (U~ 1 S~ 1 U~ 1 )\y, = U^d^S^l^U -1 , as the term localized 
on E is only generated when <9 M acts on S^ 1 . The gauge field = A^(A, S) 



in (11) enjoys the following properties, 

A" {A, S) = A(A, US) , A(A, S) = A(A^, SU' 1 ) . (14) 

Then, in terms of the A, S variables we have a double redundancy, the 
usual one associated with invariance of the Yang-Mills action under gauge 
transformations, A^ ■ T = UA^ ■ TU~ l + ^Ud^U" 1 , represented by S — > US, 
and other originated from the different ways to express the same vector field, 
combining the transformation A^ ■ T = UA^ ■ TU^ 1 + -Udjj~ x , together 
with a right multiplication of S. 

At this point, we would like to emphasize that a nonperturbative defi- 
nition of the path integral in Yang-Mills theory is still lacking. This comes 
about as a gauge fixing procedure generally leads to Gribov copies [33] in that 
regime, so that it is difficult to define an appropriate object where each phys- 
ical situation is counted only once. The restriction to the modular region has 
been usually implemented by means of the Zwanzinger action |34| . In this 
framework, in the infrared regime, the path integral has been shown to be 
dominated by configurations near the Gribov horizon. On the other hand, as 
is well-known, configurations containing magnetic objects proliferate at the 
horizon |35j-|37j. From this perspective, it is natural to fix the redundancy 
by introducing the identity 1 = A_pp[^4.] J[dU] 5[f(A u )}, in the perturbative 
sector where the Faddeev-Popov procedure is well defined. In addition, as the 
S sector parametrizes correlated monopoles and center vortices, it represents 
configurations at the horizon, relevant to describe the large distance physics. 
Giving a configuration S, gauge fixing amounts to choose a representative of 
the orbit U S. Any condition imposed on A(A, S) will be invariant under the 



[/-transformations in eq. (14). This is also the case for conditions depending 



on I(S), as it is invariant under right multiplication (cf. eq. (13)). 

In this article we shall not attempt to derive a precise construction of 
the integration measure. Rather, having the previous remarks, notation, and 
conventions in mind, we argue that it is quite natural to propose the following 
path integral, 



J YM 



j [dA][dS] A FP [A]5[f(A)]e- s ^ , (15) 
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where dS, represents the ensemble integration over monopoles and thin cen- 
ter vortices, that is supposed to include its own appropriate gauge fixing 
condition. Note that, in the trivial sector, where S = S r is regular, we have 



A(A, S r ) = A Sr and the associated contribution to (15) is the usual, pertur- 
bative one. Only in that sector A(A U , S) = A(A, SU) may be identified with 
a gauge transformation. 

As a final step, and as a guide to the construction of the effective model, 
we rewrite the partition function in the equivalent form: 

Z YM = [ [dA][dS][d\] AM^S\f(A)]e-^^ x ^^>-^^], (16) 



where = le^pj-"" an d we have introduced a color-valued auxiliary field 
A" to deal with a first-order version of 



3 The effective theory 

Let us now derive a non Abelian effective field theory for the sector of de- 
fects. The derivation will become possible by relying on the symmetries 
exhibited by the ensemble integration. This effective theory shall contain 
mass parameters, which we assume are originated from those present in a 
(phenomenological) ansatz for the action of the defects. In this regard, we 
note that up to now we have considered thin center vortices, parametrized as 
in Q. However, lattice simulations [18] point to the idea that they become 
thick objects, characterized by some finite radius of the order of lfm. More- 
over, as discussed in [26j, the stable objects in the continuum could in fact 
correspond to some deformation of the thin objects given in where the 
"thin" quantities CfJ are replaced by some smooth finite radius profiles C^. 
If this is assumed to be the case, rather than eqs. 0, g, the Yang-Mills 
action would have the form, 

S YM = J d 3 x \(^JA) ~ ^(C)) 2 + K , (17) 

where 1Z vanishes for thin center vortices. Note that the first term can be 
linearized, as we did before, by introducing the fields \ a ^. Besides, at large 
distances, approximating C° by C", this term shall originate the terms ap- 



pearing in the exponent of eq. (16), when the center vortices were considered 
to be thin. On the other hand, the second term (TZ), will be concentrated 
on the center vortices and at large distances will produce instead an addi- 
tional action Sd for the defects. Therefore, in the general case, the ensemble 
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integration must be written in the form, 



e 



Sv, m [X] = h d (f\ e -S d +if<Px\*F*(C) 



The second term in the exponent above has a local 50(3) symmetry under 
right multiplication: 5 — > SU, changing the color basis from h a -T = 5T a 5 _1 



to h' a ■ T = SUT a U~ l S~ l , that is, h' a = R(S)R{U)e a . Note that, using (|lO|), 
and that R(U) contains no defects, we have, 

7£,{C) = ^{R{U)M a R-\il)R-\S)[d,,d v }R{S)). (19) 

In other words, a regular local rotation of A™ can be translated to a regular 
local transformation of 5. Then, if 5^ were nullified, that is, if we were 
dealing with strictly thin center vortices, 5 v m [A] would be invariant under 
local 50(3) rotations, as the transformation 5 — > SU could be absorbed by 
the integration measure dS. In this regard, we would like to underline that 
this measure is to be accompanied by an appropriate gauge fixing condition 
that is invariant under right multiplication (see the discussion at the end of 
the previous section). However, in 5„ jm [A], that symmetry will be broken to 
a global one because of the thick character expected for center vortices. To 
have a simple picture about this statement we note that an action for thick 
center vortices will typically contain a Nambu-Goto term plus other terms 
describing the center vortex rigidity \27\ 138 j . These pieces can be generated, 
for instance, from a large distance approximation of the more symmetric term 
(in color space), 



J d 3 xd 3 y^(C)\ x G M (x-y)^(C)\ 



v 



where Gm is a kernel localized on a scale 1/M. Now, as we have seen in 
eq. (19), the field strength J-^ V {C) will rotate under local U transforma- 
tions. Therefore, as for any finite M the integrand above depends on the 
field strength at different spacetime points, it will change under the local 
transformations, only leaving a symmetry under the global ones. 

Based on purely geometrical/mathematical grounds, the possible kinds 
of defects can be straightforwardly classified as follows: 

i) Closed center vortices. 

ii) Monopoles and antimonopoles, joined by center vortices (each pointlike 
object is joined by a pair of center vortices). 
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iii) A particular limit of ii): A coincident pair of center vortices, which 
should correspond to an unobservable Dirac string. 

To proceed, let us consider a type ii) configuration (correlated monopoles 
and center vortices). To that end we recall that, in refs. |8]-[Tn], we have 
considered a particular case of that situation, namely, when the center vor- 
tex color points along the (locally) diagonal direction n 3 . In that case, the 
effective field describing these objects corresponds to a complex vortex field 
V. In particular, in ref. |1U| . we have shown how the ensemble integration 
over open center vortices, whose endpoints are joined in pairs to form closed 
chains, leads to an Abelian Z(2) effective theory that can be written in terms 
of V, thus making contact between the initial representation and the final 
effective field theory. For this aim, we applied recent polymer techniques 
[3"9"| HO] to deal with the end-to-end probability associated with center vor- 
tices interacting with a general vector field A M , and a scalar field needed to 
represent vortex-vortex interactions. However, it is far from straightforward 
to extend this type of derivation to the non Abelian context. Therefore, in 
our case, we will propose a model relying on the symmetries displayed by 
the initial representation, that strongly constraints the possible associated 
effective theories. 

In our case, the candidate for a vortex field has to be a real 3-component 
field (p a (a = 1,2,3), because of the global SO (3) symmetry of the action 
S[X\. We shall also introduce an isospin-2 field Q, where Q is a traceless 
symmetric 3x3 real matrix, encoding information about how the monopole- 
like defects that center vortices can concatenate are distributed. We may 
then consider in the effective theory, an invariant term Vj that couples the 
monopole and vortex sectors: 

Vi = C 4> T Q , C = constant (20) 

which is invariant under the local SO (3) transformations: 

(f)(x) ->■ R(x) <f>{x) , Q(x) -> R{x)Q{x)R T (x) . (21) 

There are also invariant terms involving just either the vortex or the 
monopole field. Regarding the former, we may include a 'potential' term V^, 
with the general structure: 

V* = ^0 T 0+^(0 T 0) 2 , (22) 

where /x and A are arbitrary constants. On the other hand, for the case of the 
monopole field, we recall that an order parameter Q, with a similar structure, 
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is well-known in the context of liquid crystals. Thus, we expect the relevant 
terms in the effective theory to be of the same kind, namely, we may include 
a potential Vq 



V Q = ^ + f A + ^ + °SA + | A 2 + FS 3 , (23) 

where A,...,F are constants, and we have introduced two independent 
SO (3) invariants ^ that can be built in terms of Q: 

5 = TrQ 2 , A = Tr Q 3 . (24) 



Thus, the three terms V^, Vq and Vj have the local symmetry (21). This 
local symmetry will be broken to its global counterpart by the kinetic terms; 
however, these terms shall be constructed in such a way that they are compat- 
ible with a local discrete gauge symmetry. This symmetry must be present, 
at least in a phase where the vacuum is symmetric (no spontaneous symmetry 
breaking) . 

In this regard, the field strength tensor JFHC) can be written as, 

— £fiupd u Cp + -e^ up e abd C^Cp, (25) 

where, 

C 1 = --h 2 -d u n 3 
9 

C 2 = — n 3 • d v fii 
9 

Cl = --h x -d v h 2 . (26) 
9 

We will show that J-^(C) can be rewritten as, 

jF;{C) =h%- h% (27) 

h a ll = ei lvp d v Cf , hl = ---e llvp h a -{d v h a xd p h a ), (28) 

Zg 

where, in the second tensor, no summation over a is understood. 



2 Being a traceless real symmetric matrix, the invariant content of Q can be generated 
by two real invariants. For example, two of its eigenvalues. 
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Let us take, for example, the third component of the field strength tensor, 

(30) 



• r it '> 'V' 



h A — p f) C A 



^3 _ „, 



-gCftvp C V C 



In order to show that 



1 



h* = -—e^ph 6 ■ {d v h 6 x dph 6 ), 
^9 



(31) 



we can simply note that 
d v h 3 = 



{fix ■ dvhzjh-i + (h 2 • d v h 3 )h 2 + (h 3 ■ d u h 3 )h 3 
(ni • d v h 3 )hx + {h 2 ■ d v h 3 )h 2 
giClhx - Cln 2 ). 



(32) 



Then, replacing in the second member of (31), and using h\ x h 2 = h 3 , etc. 



it is straightforward to make contact with (30). 

The important point is that (27) and (|28[) imply that for a fixed monopole 
background correlated with center vortices, the integral of each component 
over a closed surface dd (given as the border of a three- volume 



dS,K{C) 



1 



I 



dSp(hl-h*) = — <b dSp6p V ph a -(d u h a xdph a ), (33) 
ae ^9 Jae 



gives the II2 topological charge for the mapping dd — » h a . More precisely, 



dSp^(C) 



4ir 
9 



(34) 



where (n_(i?)) is the number of monopole (antimonopole) defects in- 

side for the component ha- 
lt may appear that the previous expression sets a preferred direction in 
color space. This impression can be dispelled by considering the effect that 



a space independent change of color basis has on the expression (34). To 
that end, we consider a new basis (n4) 3 _ i; related to the original one by a 
(constant) matrix R', namely: h' a = R'h a , a = 1,2, 3. In components, and 
using an obvious notation, the last relation means: 



(n' ) 6 = {R')bc{h a )c ■ 



(35) 



Thus, we see that: 

dSp, e^p n a ■ (d v h' a x d p h' a ) 



dSp e^p e hlh2h3 {h' a ) hl d v (n a ) h2 dp(n a ) h3 
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dSfx £(u/p ^feib2^3 (-^ )&ici (-R )&2C2 (-^ )^3C3 (^a)ci d v {fla)c2^pif^a)cz 

= (det-R') ® dSpe^ vp h a ■ (d v h a x <9 p n a ) , (36) 
where we have used the property: 

(37) 

A similar relation holds if one changes the original canonical basis by a con- 
stant rotation matrix. What this proves is that it is possible to generate a 
monopole charge along any color direction as long as one needs how to do 
that for, say, the third one. 

Thus, coming back to the discussion on the possible form of the kinetic 
terms, they must be -at least in the symmetric phase- compatible with the 
local discrete gauge symmetry: 

A« -> A£ + dpu a , (38) 

where uj a is a discontinuous function taking values ±| inside a three- volume 
■d a , and zero outside. Note also that in a phase only containing closed center 
vortices, a larger symmetry, 

K^K + V» ( 39 ) 

for any smooth (f a , is expected: in this case, the absence of monopoles would 
imply d^(C) = 0. 

Then, the kinetic terms must have the global SO (3) symmetry plus the 
local Abelian one. The simplest choice, which, in the effective theory ap- 
proach spirit we shall consistently adopt, is to minimally couple <f> and Q 

— * — * 

to A^ (note that these couplings do not have the local X(x) — > R(x)X(x) 
symmetry). Thus, the structure of the kinetic term K is as follows: 

K = K^ + Kq, (40) 

where: 

Kg = \{VpQY\VpQ) ab (41) 

where V p denotes the covariant derivative operator (consistent with the sym- 
metries mentioned above), which shall adopt a different expression when 
acting on each one of the fields. Explicitly: 

{VpQ) ab = dpQ ab -ig Q X^ acd Q dh + ig Q Q ad e dch Xl, (42) 
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where and are constants. 

The global 5*0(3) symmetry is evident, while by imposing = gg, the 
effective action will display a non Abelian gauge symmetry, and the different 
phases for the ensemble of monopoles and center vortices will arise as different 
possible vacua when the system undergoes SSB. Note also that as the field <fi 
represents center vortices that in the case of Abelian configurations posses a 
magnetic charge 2ir/g, the natural choice is g<p = 2ir/g, which also matches 



the correct dimensions in eq. (42) as [A] = 3/2, [g] = 1/2. 



Then, joining the different pieces and taking into account eq. (16), the 
following model, encoding a general ensemble of magnetic defects, can be 
proposed, 

£ eff = C v , m + ^A« + iXfflA), (43) 

C v , m = K Q + K <t> + V <j) + V Q + Vj. (44) 

We would like to underline that according to the discussion at the end of £|2j 
and be ginn ing of £j3j the symmetry displayed by the second and third terms 



in eq. (43), namely a transformation A u , accompanied by the local 5*0(3) 
rotation of A", is not the gauge symmetry that operates on A^. Therefore, 
the noninvariance of C VjTn under local 50(3) rotations of A^ is not an explicit 
breaking of the gauge symmetry in our effective model. Only in the trivial 
sector A may be associated with a gauge transformation, in other words, 
our model refers to the interaction of effective fields, parametrizing a general 
ensemble, with effective gluons represented by A^. 



4 Phase structure 

In order to analyze the possible phases of the model, it is necessary to study 
all the possible scenarios regarding both the (f) and Q dependent potentials, 
Vtj, and Vq, as well as the interaction Vj. This will yield information about 
the possible translation invariant configurations that will determine the prop- 
erties of each phase. Non translation invariant configurations, on the other 
hand, are important to understand the mechanism driving the phase transi- 
tions between them. Of course, that will require the inclusion of the derivative 
terms into the game. 

Thus we consider the minima of 

V T = V^ + Vq + Vj. (45) 

This analysis is greatly simplified if we note that Q can always be diagonalized 
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by a similarity transformation Q = R T DR, with 



-i - 2 

£> = | 'o 2 -f + 2 | . (16) 


Defining R(j) = if), the potential Vr adopts the form, 
Vt = ^+|a+^ 2 + ^A+|a 2 + ^^+^(^) 2 +C^^, (47) 



5 = (3g 2 + r? 2 )/2 
A = 3g(g 2 -?? 2 )/4 

VDiP = -|(A 2 + ^ 2 ) + |(^ 2 -^) + #3 2 - (4? 



Here, the term <5 3 that was present in eq. (23) has been discarded, as it does 
not modify the qualitative structure of the minima |32| . 

Now, to find the minima of the potential, we will suppose that the chain 
of spontaneous breaking of the symmetries is dominated by the monopole 
sector. Concretely, this approximation amounts to finding the minima of Vq, 
and using the configurations go, "go that Q adopts in those minima as a fixed 
background where we look for the vortex field configuration that minimizes 
the remaining potential. Then, the whole space of minima is generated by 
means of i?-rotations of the former. 

The minima of Vq are determined by: 

d q V Q \ =0 , SLVqI =0 , (49) 

plus the usual conditions on the second derivatives. We will consider CE > 
6D 2 /25, and will follow the discussion in |32| . where the different kinds of 
minima are obtained by varying A and B. Changing the independent vari- 
ables q and r], the region 5 3 > 6A 2 is mapped, and the strict inequality 
occurs for rj ^ 0. Then, the points obtained by simply minimizing with re- 
spect to 5, A as independent variables (in this case the potential contains a 
positive definite quadratic form) can only correspond to r] ^ 0. Otherwise, 
the potential must be minimized with the constraint <5 3 = 6A 2 , in which case 
two different situations are obtained, go — 0, i]o = or go 7^ 0, 770 = (when 
CE < 6D 2 /25, only the two latter possibilities can be realized). 

Then, the ip-&e\d minima follow from the study of the 'effective potential' 
V (?/>), defined by: 

v(V0 = ^ + VK<7o,wV), ( 5 °) 
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Figure 1: A-B phase diagram for the monopole sector, when CE > 6D 2 /25, 
D < 0. 

and which explicit form is: 

Ira., v ,, 1 



+ i0* 2 + 2C(R,)^ + ^i a +^ + ^f) 2 . (51) 

It is now clear what kind of vacua may emerge, depending on the relative 
values of the parameters. We first note that stability requires A > 0. For 
D < 0, the monopole phase diagram is that of fig. [T] [32j. If the parameters 
A, B are initially in region I, we have Q = (qo = r]o = 0), and the effective 
vortex potential results, V 1 ^) = + \ {' l P 2 ) 2 - Then, if fj 2 > 0, the mini- 
mization gives if} = 0. With this vacuum, S VtTn displays a non Abelian gauge 



symmetry, much larger than the Abelian symmetries in eqs. (38), (39), typ- 
ically obtained when monopoles and center vortices are present. Therefore, 
this phase represents a situation where monopoles and center vortices do not 
proliferate (deconfining phase). Still in the Q = phase, but with fi 2 < 0, 
the system undergoes SSB leaving an Abelian symmetry. If the mass scale 
generated for the off-diagonal fields is large, they will be suppressed and 
then the effective theory will essentially be invariant under Abelian gauge 
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transformations of the form, 



A M ■ 0o ->■ -V ■ 0o + <9 M <£>. 
Thus, recalling eq. (39), this phase describes an ensemble of closed center 



vortices. 

Now, in order to continue the analysis, it is convenient to define a complex 
field V = ttjC^i +^2), and rewrite eq. (51 ) in the form (we consider £ < 0), 

V(0) = (^ + |Cko)^ + i|C|r / o(K 2 + y 2 ) 

+ ^-2\C\q )^ + \(vV + ^i] . (52) 

When A is lowered from positive to negative values, after a first order tran- 
sition, we will enter the uniaxial nematic phase iV + (go > 0) or the N~ 
(go < 0), depending on whether B < or B > 0. These phases are charac- 
terized by t]q = 0. In what follows, to simplify the analysis, we will suppose 
fi 2 > 0. Then, if we enter the N~ phase, after a discontinuous transition, 
the effective potential V~(ip) will be minimized by ^3 = 0. In the monopole 
sector, the vacuum will be invariant under rotations around the third axis, 
while in the l/-sector this symmetry will undergo a U(l) SSB or not, de- 
pending on the sign of (/i 2 + |C|<7o)- 111 addition, the N~ phase will induce 
a mass of order g 2 , for the charged dual vector fields A* and A 2 , originated 



from the covariant derivative of Q in eq. (42). If we assume this mass to be 
large when compared with the other mass scales in the problem, these dual 
vector fields will become suppressed. 

If we further diminish A, after a second order phase transition, we will 
eventually reach the biaxial phase iV b where r] 7^ 0. As this transition is 
continuous, and we are approaching from the N~ phase, we will start with 
■03 and T] small. In the Nj, phase, the U(l) symmetry of the effective action in 
the former iV~ phase will be broken to a discrete one under 7r-rotations along 
the third axis. Again, in the monopole sector the vacuum is invariant, while 
in the vortex sector it will display SSB of the discrete 7r-rotations depending 
on the sign of (fi 2 + |C|<?o)- When this quantity is negative, at the minima, 
the V field can take a pair of values Vq, — Vq connected by a Z(2) symmetry. 
That is, the obtained effective potential coincides with the confining phase 
of the vortex model introduced by t' Hooft, relying on the possible nontrivial 
vortex correlators in the initial theory. In this phase, the spontaneous Z(2) 
symmetry breaking leads to domain walls attached to Wilson loops, thus 
providing an area law. Still in the (fi 2 + ICI^o) < case, in the intermediate 
N~ phase, the vacuum no longer displays the Abelian symmetry present 
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in the initial phase, where center vortices are only closed objects, nor the 
discrete symmetry of the last phase, typical of open center vortices whose 
endpoints are joined in pairs to monopole-like objects that proliferate. From 
this perspective, we speculate that the N~ phase might be associated with 
one where monopoles and antimonopoles are still bound in pairs. 

5 Conclusions 

We have constructed a novel non Abelian effective model for 577 (2) QCD in 
Euclidean three-dimensional spacetime that allows for the description of a 
phase diagram with a rich structure. The construction is based on a special 
parametrization of the gauge field configurations in terms of a vector 
field representing a topologically trivial sector of smooth fluctuations, 
and a local color frame n a containing defects, the nontrivial sector describing 
monopoles and thin center vortices. The frame can be written as a local 
5*0(3) rotation R of the canonical basis e a , which can be also expressed in 
the form R = R(S), where 5 is in the fundamental representation. 

This parametrization is used to write the Yang-Mills action, what defines 
the weight assigned to each configuration. On the other hand, as in any 
non-perturbative definition of the functional integration measure in a non 
Abelian gauge theory, one is faced with the usual stumbling blocks, related to 
the Gribov problem. We do not attempt to tackle this problem; rather, since 
we use the functional integral just as a guide for the subsequent derivation 
of the effective model, we use instead a definition of the measure which: (a) 
reduces to the proper one for topologically trivial configurations and (b) is 
consistent with (although not uniquely determined by) the properties of the 
gauge field parametrization used. 

The next step in the construction of the effective model proceeds with the 

— * 

introduction of an auxiliary field A M that linearizes the Yang-Mills action, and 
the incorporation of a phenomenological weight Sd that senses the geometry 
of the defects. It is at this point where the real reduction to an effective 
theory is implemented. Indeed, the symmetries are identified here, for a 
given classification of defects, what allows us to construct an effective model. 
If S d were nullified (thin objects), the partition function for the sector of 

— # 

defects should be invariant under local SO (3) rotations of A, as they could be 
absorved by a frame redefinition, transforming S under right multiplication 
by an appropriate regular 577(2) matrix U" 1 . In the Yang-Mills partition 
function, the symmetry should also be accompanied by the transformation 
A^. However, this symmetry is the one associated with the many different 

ways a given gauge field containing thin defects can be decomposed, so 
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that it is expected to be broken as soon as center vortices become thick. 
Alternatively, this could be seen as the noninvariance of the effective phe- 
nomenological action under local frame rotations, only leaving a global 
SO(3). 

An interesting point is that in order to guide the construction of the 
effective model for the ensemble integration, not only the global 5*0(3) sym- 
metry is important but also a new symmetry comes into play. At least in the 
symmetric phase, due to the topological structure of monopoles, the model 
should be invariant under a local discrete gauge symmetry. This led us to 
propose a non Abelian model describing the interaction of the natural or- 
der parameters for monopoles and center vortices with the effective gluon 
field A^. As center vortices can be attached in pairs to the non Abelian 
monopoles, the corresponding order parameters are given by fields and Q, 
carrying isospin one and two, respectively. The effective character of the 
gluons is due to the fact that gauge transformations of the Yang-Mills fields 

— * — * 

act as a left multiplication of the S sector, leaving A^ invariant. 

The effective model we introduced exhibits a rich phase diagram. For 
instance, the monopole sector of the effective potential depends on two in- 
variants, 5 = TrQ 2 , A = TrQ 3 . If this sector is supposed to dominate 
the transitions, the phase diagram inherits, by construction, some of the 
properties found in liquid crystals. In this case, if the quadratic form in 
the quantities 5 and A is positive definite, and the coefficient of the linear 
A-term is positive, an interesting chain of phase transitions is obtained. 

Initially, when the coefficient of the linear 5-term (A) is varied from posi- 
tive to negative values, a first order transition from the isotropic deconfining 
phase to a uniaxial monopole condensate takes place. In this process, in the 
vortex sector, the "third" component becomes suppressed, while the other two 
components can be arranged as an Abelian complex vortex field V display- 
ing U(l) SSB. In this example, the vortex mass scales have been supposed to 
be negligible when compared with those generated in the monopole sector. 
The further reduction of A produces a second order phase transition, and 
the monopole condensate becomes biaxial. Here, center vortices are left in a 
global Z(2) SSB phase, thus making contact with the 't Hooft vortex model, 
and arriving to the confining phase expected in 3D Yang-Mills theories. 
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